Electric field induced inversion of the sign of half-integer disclinations in 2D nematic liquid crystals 
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We study the effect of the rotation of an external electric field on the dynamics of half-integer disclination 
networks in two dimensional nematic liquid crystals with a negative dielectric anisotropy using LICRA, a Liquid 
CRystal Algorithm developed by the authors. We show that a rotation of n of the electric field around an axis of 
the liquid crystal plane continuously transforms all half-integer disclinations of the network into disclinations 
of opposite sign via twist disclinations. We also determine the evolution of the characteristic length scale, thus 
quantifying the impact of the external electric field on the coarsening of the defect network. 



I. INTRODUCTION 



II. THE MODEL 



o 

> 
O 



O 

a 

a 

o 
o 



> 

00 

in 
in 



X 



Topological defects play a fundamental role in condensed 
matter physics dEl and cosmology [3 ]. Liquid crystals are 
an example of a very rich environment where the dynamics 
of topological defects can be realized experimentally at rel- 
atively low costs, making them ideal laboratories for testing 
different scenarios for defect formation and evolution l4lfT2ll. 
External fields, such as electric and magnetic fields, can have a 
significant impact on the orientational order of liquid crystals 
(13HI11 thus affecting the coarsening dynamics of disclination 
networks. Recently it was shown that the orientation of the ap- 
plied electric field and the sign of the LC dielectric anisotropy 
may be used to control the type and topological charge of 
disclination networks fT8ll . In ifTTl a continuous transforma- 
tion of some —1/2 wedge disclination lines into 1/2 ones was 
numerically simulated on a cholesteric blue phase of a chiral 
liquid crystal, through the application of a constant electric 
field oblique to the line wedge. This provided a numerical re- 
alization of a continuous transformation between half-integer 
disclinations with topological charge of opposite sign 0. 

In the present work we provide a numerical realization of 
the inversion of the sign of all half-integer disclinations on a 
2D nematic LC with a negative dielectric constant induced by 
the rotation of an external electric field. In our implementa- 
tion all half-integer disclinations of the network transform into 
disclinations of opposite sign via twist disclinations. We also 
consider other dynamical effects associated with the applica- 
tion of the extern electric field on the nematic, by comparing 
the coarsening of the network under the rotation of the electric 
field with its evolution under a transformation of the director 
profile implemented by hand. 

The paper is organized as follows. In Sec. [TTJ the equations 
describing the relaxational dynamics of nematic LC, in terms 
of a symmetric, traceless order parameter Q a p, are presented. 
The numerical techniques are discussed in Sec. Ill In Sec. 



[TV] the effect of the rotation of an external electric field on 
the evolution of half-integer disclination networks on a two- 
dimensional nematic LC with a negative dielectric anisotropy 
is studied in detail. The conclusions and final remarks are 
presented in Sec. [V] 



The orientational order of a nematic LC without intrinsic 
biaxiality is described by a symmetric traceless tensor, Q a p, 
at every point in space, whose components are given by Q 

3 11 

Qa/3 = -S{n a n(3 - -5 a p) + -T (I dp m a mp) , (1) 

where the unit vector n is the director, determining the local 
average orientation of the molecules, 1 is the codirector, asso- 
ciated with the direction of orientational order perpendicular 
to n and m = n x 1. The variables S and T represent the 
strength of uniaxial and biaxial ordering, respectively. The 
values of S and T may be found by the diagonalization of the 
matrix 



Qott 






-(S-T)/2 
S 



(2) 



in a coordinate system where n = (0, 0, 1), 1 = (0, 1, 0) and 
m= (1,0,0). 

Static equilibrium can only be reached for a minimum value 
of the free energy (SF/SQ a p = 0). However, the evolution 
of the order parameter Q a p from a given set of initial condi- 
tions is not fully specified by the free energy functionals, and 
further assumptions have to be made on how the minimization 
process will take place. In the absence of thermal fluctuations 
and hydrodynamic flow, the time evolution of the order pa- 
rameter is given by lfT9l 



Qap(r,i) 



5F 



(3) 
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Here the dot represents derivative with respect to the physical 
time, t, and the tensor 



(4) 



^olPHv — r ( S afjL 5p u + 5av5/3fi ~ -^S a pS^ u 



satisfies T n 



and T n 



thus en- 
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suring that the order parameter Q a p remains symmetric and 
traceless. In the following we shall assume that the kinetic 
coefficient, T, is a constant. 
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Figure 1: (Color online) Projection of the director n onto the two 
dimensional grid along with the value of the parameter S at every 
grid point, at t = 499. Only half-integer disclinations remain due to 
application of an external electric field along the z direction. 



The free energy can be written as 

F = J d 3 r(T b + T el +T E ) . 



(5) 



The first term, T b , is the bulk free energy density. It describes 
the nematic-isotropic phase transition and it may be obtained 
form a local expansion in rotationally invariant powers of the 
order parameter 



TrQ 2 
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(6) 



The second term, T e \ , is the elastic free energy density. Using 
the one elastic constant approximation the elastic free energy 
density is given by 



R 



(7) 



where the constant R is the single elastic constant. The last 
term in the free energy functional, Te, is the contribution of 
the effect of an external electric field E, 



te 



(8) 



where te = 2Ae and Ae is the dielectric anisotropy. The 
molecules of a LC with positive dielectric constant, Ae > 0, 
tend to orient parallel to the external electric field, E. On 
the other hand, if the LC has a negative dielectric constant, 
Ae < 0, then the molecules tend to be align perpendicularly 
toE. 

With this free energy the equation of motion can be written 
more explicitly as 



Qap= -T[(A + CTrQ 2 )Q a p + B Q ai Q lP 



(9) 



Figure 2: (Color online) Projection of the director n onto the two 
dimensional grid along with the value of the parameter S at every 
grid point for a single snapshot of simulation B at t = 650, after the 
rotation of the electric field by an angle of tt/2. Only half-integer 
twist disclinations appear in this snapshot. 



where a comma denotes a partial derivative and X a p set- 
ting an arbitrary real matrix, X a p, to be a symmetric traceless 
matrix by the following procedure 



77 * 



(10) 



III. NUMERICAL IMPLEMENTATION 

All the numerical simulations were performed with LICRA 
(Liquid CRystal Algorithm), a publicly available set of C 
codes and MATLAB/OCTAVE routines used to solve Eq. (|9j 
using a standard second-order finite difference algorithm for 
the spatial derivatives and a second order Runge-Kutta method 
for the time integration. This software is free and it is avail- 
able at http://faraday.fc.up.pt/licra |[T8l . 

The order parameter Q has 5 degrees of freedom associ- 
ated with S, T, n and 1 (n accounts for two degrees of free- 
dom). The initial conditions for S and T are randomly gen- 
erated, at every grid point, from uniform distributions in the 
intervals [0, 2/3] and [0, 35], respectively. The director, n, is 
also randomly generated, at every grid point, using the spher- 
ical vector distribuitions routines in the GNU Scientific Li- 
brary (GSL) and the codirector, 1 , was calculated by randomly 
choosing a direction perpendicular to n. The eigenvalues and 
eigenvectors in LICRA are computed using the library GSL. 

In order to determine the evolution of the characteristic 
scale of the network LICRA numerically calculates the cor- 
relation function in Fourier space 



P(k,t) 



QapQ£,t)Qp a (-k,t) 
j rf 3 kQ^(k,t)Q^(-k,t) 



(11) 



after setting the infinite wavelength mode Q a p(0) to zero. 
Here Q a p{k) is the Fourier transform of Q a p{r) (k is the 
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Figure 3: (Color online) Similar to Figs. [2]and[T]but now at t = 900, 
(t > toff). The sign of all topological charges were inverted with 
respect to Fig. [T] 

wavenumber) and was calculated by using the library Fastest 
Fourier Transform in the West. The characteristic scale, L, 
can then be defined as 

^ = <^ 2 > = E k2p (^ *)/ E p ( k > *) • < 12 ) 

k k 

In this paper we have used LICRA to perform high res- 
olution numerical simulations of the dynamics of a texture 
network in a uniaxial nematic two-dimensional LC under an 
external electric field. All simulations were performed on a 
2048 2 grid (in the xy plane) with parameters Ax = Ay = 1, 
At = 0.1, A = -0.1, B = -0.5, C = 2.67, R = 1.0, 
r = 1.0, |E| = 0.2/v^Ae] and Ae = -0.04. 

IV. RESULTS 

The evolution of Q a p starts with the external electric field 
switched off. An electric field perpendicular to the plane of 
the simulation (along the z direction) is connected at t = 300 
and disconnected at t = 350. This procedure ensures that 
only half-integer disclinations remain in the simulation ifTSl . 
After that, three different simulations have been performed: 
simulation A with no external electric field, simulation B 
where a rotating external electric field has been applied be- 
tween t on = 500 and t Q ff = 800 and simulation C where a 
modification by hand of the director and co-director profiles 
was made at t = 500. Fig. [T] shows the value of the order 
parameter S and the projection of the director n onto the two 
dimensional grid for t = 499, which is equivalent for simula- 
tion A, B and C. 

In simulation B the components of the electric field are 
given by 

E x = 0, 

E y = |E|sin(<p), 

E z = |E|cos(y>), (13) 



Figure 4: (Color online) Projection of the director n onto the two 
dimensional grid along with the value of the parameter S at every 
grid point of a single snapshot of simulation C at t = 900, after 
changing the director configuration profile by hand according to Eq. 
|T6| . The sign of all topological charges were inverted with respect 
to Fig. [j] 

with 

^=3^-500). (14) 

The electric field was applied at t on = 500 along the z di- 
rection and then rotates around the x direction by an angle of 
cp, on the yz plane. At each time step (At = 0.1) the angle, 
ip, was increased by 7r/3000 radians until a total angle of tt 
(at t ff = 800). The negative dielectric anisotropy, Ae < 0, 
ensures that the director remains perpendicular to the field en- 
suring that n • E = 0. If the defect centers were static then, at 
each time- step, the components of the director would be equal 
to 

Tlx T^xo 1 

n y = n^cosO), 

n z = -n yo sm(<p), (15) 

where (n XQ , n yo , 0) are the components of the director at 
t = 500 in simulation A. In practice this result must be com- 
plemented with the coarsening dynamics of the network. 

Fig. [2] presents the order parameter S and the projection of 
the director, n, onto the two dimensional grid for a snapshot 
taken from simulation B, at t = 650, when ip = 7r/2. In 
this case, the director is given by n = (n Xo ,0, —n yo ) and 
twist disclinations are formed on the xz plane. Note that the 
presence of the electric field increases the size of the defect 
cores. 

Fig. [3] shows a snapshot of simulation B for t = 900 af- 
ter the electric field has been rotated by an angle of 7r and 
then switched off when at t = 800. In this case, the compo- 
nents of the director are given by n = (n Xo , — n yo , 0) and all 
the half-integer disclinations are back into the xy plane, with 
an opposite topological charge to the one they had in Figure 
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Figure 5: (Color online) Evolution of the characteristic length, L(t), 
for simulations A, B and C. Since the curves for simulations A 
and C are almost indistinguishable we intercalate the corresponding 
datapoints in order to improve the visualization. 

[T] A similar realization of the exchange of sign of all half- 
integer disclinations network is possible through the rotation 
of the electric field around any other axis of the simulation 
plane. Besides this transformation of the topological charge, 
the electric field produces other effects on the network coars- 
ening that need to be addressed in more detail. 

In simulation C we have implemented by hand the modi- 
fication of the director and co-director profiles corresponding 
to a rotation of 7r of an external electric field. This transforma- 
tion was performed with the purpose of isolating the effect of 
the inversion of the sign of the half-integer disclinations from 
the modification of the coarsening dynamics associated with 
the presence of an external electric field. We have modified 
the components of n and 1 of the configuration at t = 500, 
shown in Fig. [T] so that 

lx — n y0 > ly ~ nxO j (16) 

where (n XQ , n yo , 0) are the components of the director n pre- 
sented in Fig. [T] The results for the projection of the director 
onto the lattice plane are presented in Fig. [4] along with the 
order parameter S for a snapshot of simulation C at t = 900. 
Fig. [4] is very similar to Fig. [T] except for the inversion of the 
sign of the topological charges and the coarsening which took 
place between t = 499 and t = 900. 

The evolution of the characteristic length, L(t), of simu- 
lations A, B and C is presented in Figure [5] Given that the 
curves for simulations A and C are almost indistinguishable 
we intercalate the corresponding datapoints in order to im- 
prove the visualization. The effect of the rotation of the elec- 
tric field is visible in simulation B, between t = 500 and 
t = 800. An acceleration of the coarsening, reflected in 
the larger slope of L 2 (t) may be observed between t = 500 
(cp = 0) and t = 650 (cp = tt/2). It is also clear that the mem- 
ory of the connection of the field is preserved by the coarsen- 
ing dynamics in simulation B, even after the electric field is 



Figure 6: (Color online) Total free energy, given by Eq. j5j, for a 
snapshot of simulation A at t — 900. The corresponding image for 
simulation C is identical. 

disconnected. This is reflected in the larger value of the char- 
acteristic length scale of the network at late time, compared to 
simulations A and C. 

The comparison of the evolution of L 2 (t) for simulations 
A and C demonstrates that modifying the director and co- 
director profiles by hand using Eq. fl6| ) does not modify the 
coarsening dynamics. This happens because the bulk and elas- 
tic components of the free energy defined by Eqs. ^ and 
([7]), are invariant by a rotation of the directors around an ar- 
bitrary axis. Such rotation originates a rotated director and 
co-director network profile with the same free energy at every 
point in space. This result can be confirmed in Fig. [6] which 
shows the total free energy, given by Eq. ([5]), for a snapshot 
of simulation A at t = 900 (the corresponding image for sim- 
ulation C is identical). As expected, the free energy profile is 
axially symmetric around the defect centers. Hence, the ac- 
celeration of the defect coarsening due to the application of 
an external electric field is associated with the modification 
of the degree of alignment of the molecules, reflected in the 
modification of the values of the order parameters S and T, 
and it is therefore insensitive to the inversion of the sign of the 
topological charges. 



V. CONCLUSION 

In this work we have numerically demonstrated that a ro- 
tation by 7r of an external electric field around an axis of the 
plane of a 2D nematic liquid crystal induces an inversion of 
the sign of all half-integer disclinations, when a negative di- 
electric anisotropy is considered. We have further analysed 
the the impact of the exernal electric field on the coarsening 
dynamics. 

The inversion of the topological charge investigated in the 
present paper may be realized experimentally in a simple 
setup. The topological arguments do not depend on the model 
parameters and apply to the different phases of a liquid crys- 
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tal. Although we have assumed that the inversion was induced support, 
by the rotation of an external electric field, a magnetic field 
induced transition is also possible in the case of negative dia- 
magnetic anisotropy. 
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